Intermediate phase in the spiral antiferromagnet Ba2CuGe207 
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The magnetic compound Ba2CuGe2 07 has recently been shown to be an essentially two-dimensional 
spiral antiferromagnet that exhibits an incommensurate-to-commensurate phase transition when a 
magnetic field applied along the c-axis exceeds a certain critical value H^. The T = Q dynamics 
is described here in terms of a continuum field theory in the form of a nonlinear a model. We are 
thus in a position to carry out a complete calculation of the low-energy magnon spectrum for any 
strength of the applied field throughout the phase transition. In particular, our spin-wave analysis 
reveals field- induced instabilities at two distinct critical fields Hi and H2 such that H\ < He < H2- 
Hence we predict the existence of an intermediate phase whose detailed nature is also studied to 
some extent in the present paper. 



I. INTRODUCTION 

A recent experimental investigation§~il of the mag- 
netic properties of Ba2CuGe207 in its low-temperature 
phase {T < Tn — 3.2 K) estabhshed the occurrence of 
spiral antiferromagnetic itppder due to a Dzyaloshinskii- 
Moriya (DM) anisotropyH^Q. A schematic illustration of 
the spiral abstracted from experiment may be found in 
Fig. 5 of Ref. 1. li was further demonstrated that 
a Dzyaloshinskii-typeO commensurate-incommensurate 
(CI) phase transition is induced by a magnetic field H 
applied along the c-axis. As the field approaches a crit- 
ical value i?c ~ 2 T, the spiral is highly distorted while 
its period (pitch) grows to infinity. For H > He the 
ground-state configuration is thought to degenerate into 
a uniform spin-flop state. This phase transition is sim- 
ilar to the cholesteric-nematic transitioii mduced by an 
external magnetic field in liquid crystalsuEj. 

It is of obvious interest to describe theoretically 
the magnan excitations measured by inelastic neutron 
scatteringEl, but progress has been hindered by the great 
formal complexity of the calculation. Here we explore 
a new approach in which the original discrete system is 
replaced by a continuum field theory. We are thus able 
for the first time to carry out a complete calculation of 
the low-energy excitation spectrum for any strength of 
the applied field and any direction of spin-wave propa- 
gation. In addition, our analysis reveals the existence of 
a new intermediate phase whose properties we examine 
and compare with experiment. 

In Sec. II the low-energy dynamics is described in 
terms of a nonlinear a model that is compatible with 
symmetry. In Sec. Ill we present a brief demonstration 
of the conventional CI transition which will provide the 
basis for all subsequent work. The complete field theory 
is first applied in Sec. IV for an analytical calculation 
of the field dependence of the magnon spectrum in the 
high-field commensurate phase. Interestingly, the uni- 
form spin-flop state is shown to be locally stable only for 



H > H2 > He where the new critical field H2 is predicted 
to be equal to 2.9 T. A first contact with the measured 
spectrum is also made in Sec. IV. 

The main thrust of our calculation is presented in 
Sec. V where the determination of the magnon spec- 
trum in the low-field spiral phase is reduced to a quasi- 
one-dimensional band structure problem ±hat is solved 
numerically. While an earlier calculationO of the spec- 
trum at iJ = is confirmed, we are also in a position 
to analyze existing experimental data at nonzero field 
and to predict the results of possible future experiments. 
A byproduct of this analysis is yet another critical field 
Hi = 1.7 T < He beyond which the flat spiral ceases 
to be locally stable. Therefore, the combined results of 
Sees. IV and V suggest the existence of an intermediate 
phase in the field region Hi < H < H2 whose nature is 
studied in Sec. VI where we show that a nonflat spiral 
becomes energetically favorable. The main results are 
summarized in the concluding Sec. VII, while discussion 
of some technical issues is relegated to two Appendices. 

II. LOW-ENERGY DYNAMICS 

The unit cell of Ba2CuGe207 is partially illustrated in 
Fig. 1 where we display only the magnetic Cu sites. The 
lattice constants are a = b = 8.466 A and c = 5.445 A. 
Since the Cu atoms form a perfect square lattice within 
each plane, with lattice constant d ~ a/\/2 w 6 A, it 
is also useful to consider the orthogonal axes x,y and z 
obtained from the original crystal axes a, h and c by a 
45° azimuthal rotation. The complete magnetic lattice 
is formally divided into two sublattices labeled by A and 
B because the major spin interaction between nearest 
in-plane neighbors is antiferromagnetic. In contrast, the 
interaction between out-of-plane neighbors is ferromag- 
netic and weakfil. Therefore, the interlayer coupling is 
not crucial for our purposes and is thus ignored in the 
following discussion which concentrates on the 2D spin 
dynamics within each layer. 
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FIG. 1. Partial illustration of the unit cell of Ba2CuGe207 
displaying only the magnetic (Cu) sites denoted by solid cir- 
cles. 



The space group of this crystal is D^^^ or P42im and 
imposes significant restrictions on the possible types of 
spin interactions. Such -symmetry constraints underlie 
most of the earlier workETu but were not spelled out in 
sufficient detail. We have thus found it necessary to carry 
out afresh a complete symmetry analysis, including both 
nearest-neighbor (nn) and next-nearest-neighbor (nnn) 
couplings. For the moment, we restrict attention to nn 
interactions and write the 2D spin Hamiltonian as the 
sum of four terms : 



W = We + Wdm + Wa + Wz. 



(2.1) 
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FIG. 2. Illustration of the dimerization process on a finite 
portion of the 2D lattice cut along the axes x and y. The 
indices a and /3 advance along the crystal axes a and b not 
shown in this figure. The meaning of the Roman labels on 
bonds connecting nn sites is explained in the text. 



Here 



vectors (2.4) which are responsible for the observed spiral 
magnetic order or helimagnetism. 

The third term in Eq. (2.1) contains all "symmetric" 
anisotropics. Since single- ion anisotropy is not possible 
in this spin s — ^ system, the most general form of Wa 
is 

^ - ^ E E {siS'i + Si^l) ' (2-5) 

<kl> i,j 

where the indices i and j are summed over three values 
corresponding to the Cartesian components of the spin 
vectors along the axes x,y and z. Accordingly, Gu = 
{G^li) are 3x3 symmetric matrices, one for each bond 
< kl >. Again, there exist four distinct such matrices: 



--E 

<kl> 



Jkli^k ■ S;) 



(2.2) 



describes the isotropic exchange over nn in-plane bonds, 
denoted by < fc^ > , with Jki = J for all such bonds. 
Similarly, 



Wdm 



E 

<kl> 



Tfki ■ (Sk X Si) 



(2.3) 



stands for antisymmetric DM anisotropy where the vec- 
tors T>ki assume four distinct values: 



De2 
De2 



Till 



Dei+D'es, 
-Dei-D'es, (2.4) 



which are distributed over the 2D lattice as shown in Fig. 
2 where nn bonds are accordingly labeled by I, 11, III or 
IV. Here D and D' are two independent scalar constants, 
while 61,62 and 63 are unit vectors along the x,y and z 
axes of Fig. 1. It should be noted that the z-components 
of the DM vectors alternate in sign on opposite bonds, a 
feature that could lead to weak ferromagnetism. No such 
alternation occurs for the in-plane components of the DM 
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(2.6) 

which are all expressed in terms of the four scalar pa- 
rameters ifi, if 2, ifa and if 4. The latter may be further 
restricted by the trace condition ifi -I- if2 -I- ifa = be- 
cause the isotropic component of the ex change interac- 
tion is already accounted for by Eq. (E^). Finally, 



Wz = -J2{gfiBil-Si) 



(2.7) 



describes the usual Zeeman interaction with an external 
field H. 

The discrete Hamiltonian could be employed to ana- 
lyze this system by standard spin- wave techniques, but 
the calculational burden is rather significant and has so 
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far prevented a complete determination of the magnon 
spectrunJj. Nevertheless, the relevant low-energy dy- 
namics can be efficiently calculated in terms of a con- 
tinuum field theory which provides a reasonable approx- 
imation for Ba2CuGe207 because the period of the ob- 
served spiral is equal to about 37 lattice constants along 
the a;-direction. A similar approach is often. .invoked in 
the related subject of weak ferromagnetisn£3£j and can 
be implemented by a straightforward step-by-step pmpa- 
dure starting from the original discrete HamiltoniarO'Ej. 

The first step is to group spins into dimers as shown 
in Fig. 2. Each dimer contains a pair of spins denoted 
by A and B and labeled by a common set of sublattice 
indices a and (3 that advance along the crystal axes a 
and b. A more convenient set of variables is given by the 
"magnetization" m and the "staggered magnetization" 
n which are defined as 



2s 



(A + B) 



2s 



(A-B), 



(2.8) 



and satisfy the classical constraints m • n = and -t- 
= 1. We also introduce space-time variables according 

to 



Tj = V2ea, C = V2ef3, 



T^2sV2eJt, (2.9) 



where e is a dimensionless scale whose significance will 
become apparent as the discussion progresses. The final 
result will be stated in terms of the coordinates 



V2 ' 



V2 ' 



(2.10) 



along the x and y axes of Fig. 1. One should keep in 
mind that actual distances are given by xd/e and yd/e 
where d — is the lattice constant of the square 

lattice formed by the Cu atoms. Finally, we introduce 
rescaled anisotropy constants and magnetic field as 



D 

71' 



V2D' 



eJ 



2sV2ej' 



(2.11) 



where we display only those combinations of constants 
that survive in the effective low-energy dynamics. In par- 
ticular, the constant K4 does not appear to leading order. 
The further notational abbreviations 



Kq 



(2.12) 



will prove convenient in all subsequent calculations. 

Now, a consistent low-energy expansion is obtained by 
treating m as a quantity of order e while n is of order 
unity. To leading order, the classical constraints reduce 
to 



m n = 0, 



(2.13) 



m is expressed entirely in terms of n by 

m = [n X (h + d^ - n X h)] - l^in, (2.14) 
2\/ 2 2 

and the T — dynamics of the staggered magnetization 
n is governed by the Lagrangian density C = Cq — V 
where 

£0 = + h- (n X h), 

V = ^{dm - Ae2 X n)^ + ^{8211 - Aei x n)^ 



(n-h)^ + (hxd,)-n. 



(2.15) 



The dot denotes differentiation with respect to the time 
variable r, di and 82 are partial derivatives with respect 
to X and y, and (ni, 712, n^) are the Cartesian components 
of n along the axes xyz of Fig. 1. Consistency requires 
t hat a ll physical predictions derived from Eqs. (2.14) and 



( ^.15 ) must be independent of the specific choice of the 
scale parameter e. This fact will be explicitly demon- 
strated or used to advantage in the continuation of the 
paper. 

We have further examined possible modifications of 
the low-energy dynamics due to nnn spin interactions 
along the diagonals of the Cu plaquettes. Our symme- 
try analysis revealed that both antisymmetric (DM) and 
symmetric anisotropics are present over nnn bonds and 
introduce a new set of parameters. Nevertheless, in the 
continuum limit, all new parameter s me rge with those 
already present in the Lagrangian ( p^ . The implied 
remarkable rigidity of the effective low-energy spin dy- 
namics is obviously due to the special crystal structure 
of Ba2CuGe207. 

In the remainder of this section we make contact wi 
the static energy functional derived by Zheludev et al 
restricted to T = 0, which appears to differ in some re- 
spects from the potential V of Eq. ( ^.15 ). First, we 
note that we have omitted from the potential some ad- 
ditive ficld-dcpendcnt constants which play no role ex- 
cept to relate the energy to the magnetization. The lat- 
ter will be obtained in Sec. Ill by a direct application 
of Eq. (2.14). A more interesting point concerns the 
special choice of exchange anisotropy madcpin Ref. 5, 
which was suggested by the work otKaplaEllj, Shekht- 
man, Entin-Wohlman and AharonytZI, and is referred to 
as the KSEA anisotropy. If the original perturbative 
derivation of the antisymmetric DM interactionQ is car- 
ried to second ordeiO, a symmetric anisotropy re sult s 
that is described by a special case of the matrices ( |2.6| ) 
with 



K, = 0, 



^ 2J' 



D'2 
"27' 



Ka 



DD' 
~2T' 



(2.16) 
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in addition to a simple renormalizati on of the exchange 
constant J. The parameter kq of Eq. ( ^.11 ) is then given 
by kq = — X'^ and the parameter k of Eq. (2.12) van- 
ishes. Since a nonzero k is ahowed by symmetry, we shall 
keep it throughout our theoretical development. How- 
ever, our numerical demonstrations will also be restricted 
to the KSEA limit (k = 0). 

Fin ally, the term (h x d^) • n in the potential V of Eq. 
( 2jl5) is absent from the energy functional of Zheludev et 
al.EI. A contribution of that natuce is present in the early 
work of Andreev and MarchenkoE3 and plays a sigaificant 
role in various aspects of weak ferromagnetismll3. This 
term vanishes when the field is applied along the c-axis 
(h X = 0) and thus does not affect the analysis of 
the CI transition. However, such a term is important 
in the case of an in-plane magnetic field which is also 
of experimental interestEl and is briefly discussed in the 
concluding paragraph of Sec. HI. 



III. GROUND STATE 

An important first step in the calculation of the T — 
dynamics is the search for the classical spin configuration 
that minimizes the static energy 



W = / Vdxdy, 



(3.1) 



where V is the potential of Eq. ( 2.15| ). For a field applied 
along the c-axis, h = (0, 0, h), the potential is given by 



V = 



1 r 



{dmY + {d2ny + -f^nl + X^ 



-A [{diui - 82112) ns ~ {nidi - ^292) 113 ] , (3.2) 

which depends only on the parameter A that measures the 
strength of the in-plane component of the DM anisotropy, 
and the combination of parameters 

-f^ = K + X"^ + h'^ (3.3) 

that includes the external field h. A notable feature of the 



potential (3.2) is its invariance under the simultaneous 
transformations 



■ + iy ^ {x + iy) e*'^o. 



ni + in2 



{rii 



in2) 



(3.4) 



This is a peculiar realization of U{1) symmetry in that 
the usual 2D rotation of spatial coordinates with an angle 
'00 is followed by an azimuthal rotation of the staggered 
magnetization with an angle —ipo. 

The minimizaiiun problem was extensively studied in 
the earlier workEl 13. Here we briefly describe a slightly 
simplified version of the obtained solution in order to 
establish convenient notation for our subsequent dy- 
namical calculations. If we invoke the usual spherical 
parametrization of the unit vector n defined from 



rii + 177.2 = sin Q e'*. 



ris = cos 9, 



(3.5) 



the minimum of the energy is sought after in the form of 
the one-dimensional (ID) Ansatz 



e = 9{x), 



$ = 0, 



(3.6) 



which assumes that the staggered magnetization is con- 
fined in the xz-plane and depends only on the spatial co- 
ordinate X, modulo a U ( 1) transformation given by Eq. 
( ^.41 ). The potential (3.2) then simplifies to 



"4 



A) + j^cos^ 



(3.7) 



where the prime denotes differentiation with respect to 
X, and stationary points of the energy (3.1) satisfy the 
ordinary differential equation 9" + j^cos0 sinO = whose 
distinct feature is that it does not depend on A. A first 
integral of this equation is given by 9'"^ — ■y^cos'^9 = C = 
5"^, where we anticipate the fact that the minimum of 
the energy is achieved at positive integration constant 
C. Thus the desired solution = 9{x) is given by the 
implicit equation 



\/ 5'^ + 7^cos^?? ' 



(3.8) 



and is a monotonically increasing function of x. The cor- 
responding spin structure repeats itself when 9 is changed 
by an amount 27r ; i.e., when x advances by a distance 



d9 



y/S"^ + 7^cos^0' 



(3.9) 



which will be called the period of the spiral. The free 
parameter 5 is determined by the requirement that the 
average energy density w — /J^ Vdx is a minimum, 
where V is the potential (3.7) calculated for the specific 
configuration (3.8). A direct computation shows that 5 



must satisfy the algebraic equation 



d9 ^J 8^ +-1^cos^ 9 = A, 



and the corresponding energy density is 



w 



(3.10) 



(3.11) 



The configuration described above will be referred to as 
the flat spiral because the staggered magnetization is con- 
fined in the a;2;-plane. 

It is clear that the root 5 of Eq. ( 3.10| ) decreases with 
increasing 7. In fact, 5 vanishes at a critical valu e of 7 
which is easily calculated by setting (5 = in E q. ( 3.1(i| ) 
to obtain 7 = 7c = A7r/2. In view of Eq. (3_^), the 
corresponding critical field is given by 



7r2 
4 



1 A^ 



1/2 



(3.12) 
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and a spiral state is possible only for h < he- At the 
critical point, the energy density ( 3.11 ) becomes w = 
and is equal to the energy of the uniform spin- flop 
state n = (1, 0, 0). The latter is a stationary point of the 
energy functional for any strength of the applied field 
and is thought to be the absolute minimum for h > he- 
The actual stability of the spin- flop state for h > he, and 
of the spiral state for h < he, will be addressed more 
carefully in Sees. IV and V. 

Next we calculate the T = magnetizati on m — 
(mi, 7712, ms) which can be obtained from Eq. ( p. 14 ) ap- 
plied for the static configuration n = (sin0, 0, cos0) and 
averaged over the period L of the spiral. The only term 
that survives in the average is 



ma 



eh 1 



2V2L 



[ {1- cos^O) dx, 
Jo 



(3.13) 



and can be expressed in terms of quantities already con- 
sidered, namely 

For h > h e, the spin-flop state n — (1, 0, 0) is inserted in 
Eq. ( 2.14 ) to yield after a trivial computation 



ma 



eh 

2^2' 



h > he 



(3.15) 



while mi = and 1712 — —eX'/2^/2. The latter formula 
is the only place where the oscillating component of the 
DM anisotropy appears and produces a field-independent 
weak ferromagnetic moment along the y-axis. 

In order ta make definite quantitative predictions we 
use as inputEl the spin value s = 1/2, an exchange con- 
stant J — 0.96 meV, and a gyromagnetic ratio g = gc = 
2.474 for a field applied along the c-axis. Concerning 
anisotropy, we adopt the KSEA limit (k = 0) and thus 
the only relevant parameter is A which may be estimated 
from the observed spin rotation by an angle AO = 27rC 
over a distance d = a/\/2 along the x-axis. The incom- 
mensurability parameter C is related to the period L of 
Eq. (3^) by ( — e/L , w here e is the scale parame- 
ter introduced in Eq. (2.9). One may actually choose 



the free parameter e as e — D/J and thus A = 1 a nd 
= K + + h"^ = 1 + h^. At zero field, Eq. ( plo| ) is 
applied for A = 1 = 7 to yield S'^ = 0.53189772 and the 
period is calculated from Eq. (3.9) as L 
Hence, e = QL 

value C = 0.0273 measured at zero fieldnl. To summarize, 
our final choice of constants is 



6.49945169. 

0.1774, where we have also used the 



K 
„2 



0, 

l + /^^ 



A = l, (3.16) 
: = L»/J = 0.1774, 



and should be completed with the stipulation that the 
unit of field {h — 1) corresponds to 2s\/2e^/5cMs = 
1.682 T, while the unit of frequency (energy) is 2s\/2e J = 



0.5 





FIG. 3. T = G theoretical predictions for the field depen- 
dence of (a) the energy density w and (b) the incommensu- 
rability parameter C,. Solid lines correspond to the flat spiral 
constructed in Sec. Ill and dashed lines to the nonflat spiral 
calculated in Sec. VI. The three vertical dotted lines indicate 
the location of the three critical fields h\ — 1.01, he — 1.21 
and h2 — 1.73, measured in units of 1.68 T. Experimental 
data (open circles) measured at T = 2.4 K were extracted 
from Fig. 4 of Ref. 4. 

0.241 mcV . The magn etization per Cu atom is given by 
Eqs. ( 3.14 ) and ( 3.15 ) in units of sgcZ-ts = 1.237/iB. Dis- 
tance is measure d in u nits of d/e = 33.75 A. 

The constants (3.16) are inserted in Eq. ( 3.12| ) to yield 
a critical field /ic = 1.21 in rationalized units or He = 2.04 
T in physical units. This theoretical prediction is consis- 
tent with experiment and is thought to be a good in- 
dication that the KSEA limit (k ^ 0) may provi de an 
accur ate description of anisotropyO. Now, Eqs. ( |3.10| ) 
and (3.9) are applied with A = 1 and 7^ = 1 + Jrto 
yield the root S = S{h) and the period L = L{h) at field 
h. The fiel d dep endence of the energy density computed 
from Eq. ( 3.11 ) for h < he, and w — 1/2 ioi h > he, is 
depicted by a solid line in Fig. 3a. Similarly, the field de- 
pendence of the incommensurability parameter ^ — C{h) 
is calculated from 



cih) _ m 
m L{hy 



(3.17) 



where C(0) and i(0) are the zero-field parameters already 
discussed, and is depicted by a solid line in Fig. 3b. The 
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results of Fig. 3 will be completed and further discussed 
in Sec . VI. The s ame numerical data may be employed in 
Eqs. ( ^.14 ) and (3.15) to calculate the field dependence of 
the magnetization and the corresponding susceptibi lity. 



Finally, we return to the U{1) transformation 



which may be applied to the special solution (3^) to 



3.4) 



yield a family of degenerate ground-state configurations: 



Q ~ 9{x cos-00 + y auiipo), 



-Vo, (3.18) 



where ipo is an arbitrary angle. The propagation vector 
of the resulting spiral forms an angle tpQ with the a;-axis, 
while the normal to the spin plane forms an angle 7r/2— -00 
with the same axis. For the special rotation ipo = 7r/4, 
the magnetic propagation vector and the normal to the 
spin plane are parallel (screw- type spiral). This symme- 
try operation is the, basis for the bisection rule discovered 
by Zheludev et al.Q when the external field is applied in 
a direction perpendicular to the c-axis, at an angle xo 
with respect to the a;-axis. The normal to the spin plane 
rotates almost freely to align with the external field, and 
thus Xo = 7i'/2 — V'O: in order to minimiz e (eli minate) the 
positive term (n • h)^ in the potential ( 2.15 ). The new 
term (h x d^) • n in the above potential does not affect 
the bisection rule but it does modify the profile of the 
spiral. For example, when the field is applied along the 
y-axis (xo — 7''/2, ?Ao = 0, h X dz = X'hei) the staggered 
magnetiza tion is again confined in the a;z-plane but the 
potential (3.7) becomes 



9' - \f + 7'cos2 



X'h sind. 



(3.19) 



where 7^ — k + is now field independent. Never- 
theless, the external field reappears in a different form 
and requires a new calculation of the spiral based on Eq. 
(3.19). Such a calculation might actually explain the ob- 
served (weak) field dependence of the magnitude of the 
magnetic propagation vectoia and provide an estimate 
for the strength A' (or D') of the oscillating component 
of the DM anisotropy. 



IV. SPIN-FLOP PHASE 

We now begin to address questions of dynamics based 
on the complete Lagrangian C — Cq — V oi Eq. ( 2.15 ) 
applied for a field h = (0,0, fe). If we also insert the 
spherical parameters (3.5) we find that 



+ sin^e $2 ) + h sin^e 4 



(4.1) 



and 



(ve)^ + sin^e (v$)^ + -f^cos'^e + x^ 

-X [cos6 sin9 (sin$ 9i<I> -|- cos$ 82^) 
-cos$ 9ie-Fsin$ 926], (4.2) 



where V = {81,82) is the usual 2D gradient operator, 
while the Laplacian will be denoted in the following by 
l\^8l + 8l. 

We first study the high-field commensurate phase {h > 
he) where the absolute minimum of the classical energy 
is thought to be the uniform spin-flop state n = (1, 0, 0) 
or Q = tt/2 and <I> = 0. Small fluctuations around this 
state are calculated by intr oducing 8 = tt/2 + f and 
^ ~ g in Eqs. (4.1) and (4.2) and keeping terms that are 
at most quadratic in the small amplitudes / = f{x,y,T) 
and g = g{x,y,T). Linear terms do not appear because 
we are expanding around a stationary point of the en- 
ergy functional, whereas constants and total derivatives 
can be omitted because they do not contribute to the 
equations of motion. Thus the corresponding linearized 
equations are found to be 



/-A/-|-72/ = 2A52.g, 



g-l^g^-2X82f. (4.3) 



Performing the usual Fourier transformation with fre- 
quency oj and wave vector q = (91, 92) one obtains a ho- 
mogeneous system whose solution requires that the cor- 
responding determinant vanish. This condition leads to 
two branches of eigcnfrcqucncies 



'^±(q) = 



^1 



\ (7' ± y7+16A^ 



-,1/2 



(4.4) 



which will be referred to as the optical or acoustical 
mode, corresponding to the plus or minus sign, respec- 
tively. 

A notable feature of the calculated dispersions is their 
strong anisotropy. In particular, the low-q acoustical 
branch reads 



^-(q)« + -4AV7') <& 



2 /„,2\ „2 1 1/2 



(4.5) 



and demonstrates that the spin-wave velocity depends 
on the direction of propagation. It also makes it clear 
that an instability arises when 7^ < 4A^. In fact. 



the complete acoustical frequency of Eq. (4.4) becomes 
purely imaginary over a nontrivial region in q-space when 
7^ = K + A^ + /i^ < 4A^. Therefore, the uniform spin- flop 
state is unstable for h < h2 where 



h2 = \/ SA^ — K 



(4.6) 



is a new critical field. For our choice of parameters ( 3.16| ) 
/i2 = ^/3 or H2 = 1.682x/3 = 2.91 T. The important 
conclusion is that the spin-flop state is locally stable only 
for H > H2> He. 

It is also interesting to examine the gap of the optical 
branch at q = where w+(q = 0) = 7 = (^-f A^ -f ft^)^/^. 
This result may be used to illustrate our earlier claim 
concerning the role of the scale parameter e. If we recall 
the definition of the rescaled parameters ( p. 11 ) and also 
include the factor 2s^/2eJ to account for the physical unit 
of frequency, the calculated gap is independent of s and 
is expressed entirely in terms of constants that appear 
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FIG. 4. Theoretical magnon dispersions in the spin-flop 
phase, for spin-wave propagation along the a;-axis {Q2 ~ 0) 
and the j/-axis (Qi — 0). 
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in the original discrete Hamiltonian of Sec. II. Hence, in 
the KSEA hmit, we find that 



w+(q = 0) = 



-,1/2 



(4.7) 



in agreement with the magnon gap given in Rcf. 5. In- 
cidentally, this special result is the only feature of the 
spectrum actually calculated in the above reference for 
nonzero field. 

The complete dispersions are illustrated in Fig. 4 for 
_ff = 3 T, and for spin-wave propagation along the x- 
or the y-axis. The anisotropy of the spectrum is made 
especially apparent by the fact that the dispersion of the 
acoustical mode is strictly linear in the z-direction, but 
almost ferromagnetic-like in the y-direction because the 
chosen field is only slightly greater than the critical field 
H2 ~ 2.9 T. Th e nu merical data for Fig. 4 were ob- 
tained from Eq. (4.4) ap plied for our choice of units and 
constants given in Eq. ( 3. IE ). Thus we set A = 1 and 



Y = 1 + h\ with h = 3/1.682 = 1.784, and also include 
an overall factor 0.241 meV to account for the physical 
unit of energy. Finally, Q = eq is the wave vector de- 
fined on the complete square lattice formed by the Cu 
atoms within each layer, while relative units are de fined 
from Q[r.l.u.] = Q/2t: = (e/27r)g. Therefore, Eq. (HJ) 
is applied with q — (27r/e) Q[r.l.u.] — 35.418 Q[r.l.u.]. 

Unfortunately, there seem to exist no experimental 
data in thcviield region H > 3 T. In fact, the only pub- 
lished datacl were obtained for = 2.5 T < H2 and 
spin-wave propagation along the x-axis. For this spe- 
cial direction (92 = 0) the theoretical dispersions (4.4) 
do not "see" the instability. One may then deliberately 
apply them for = 2.5 T and compare the results to 
the actual data, as is done in Fig. 5 where a systematic 
disagreement is apparent in both dispersions. In particu- 
lar, the numerical fits to the data represented by dashed 
lines indicate a significant 20% reduction in the measured 
spin-wave velocity; as was already noted in Ref. 5. 

Of course, our ear lier discussion makes it clear that 
the dispersions (4.4) cannot be applied for H = 2.5 



FIG. 5. Theoretical magnon dispersions (solid lines) in 
the spin-flop phase, deliberately anplied for H = 2.5 T < H2- 
Circles denote experimental datafl and dashed lines are nu- 
merical fits to the same data. 



T because the corresponding ground state is predicted 
to be unstable. At best, the fully polarized spin-fiop 
state n = (1,0,0) survives in the field region H < H2 
as a metastable state thanks to some small tetragonal 
anisotropy that may be present in the discrete systemEl 
but drops out of the leading continuum approximation. 
An appealing scenario suggested by our calculation is 
that the system actually enters a different (intermediate) 
phase for H < H2 which consists of some sort of mixed 
domains with no definite axis of polarization. Such a 
picture could explain the effective reduction of the spin- 
wave velocity, also taking into account the anisotropy of 
the acoustical mode. 

As mentioned already, the continuum model does not 
contain anisotropics that would necessarily polarize the 
staggered magnetization along the x (or the y) axis. In- 
stead, there is a family of degenerate spin-fiop states 
n = (cos $oi sin $oi 0) with the same energy for any con- 
stant angle <i>Q . The corresponding small fiuctuations are 
now stu died by in trod ucing O — 7r/2 -f / and ^ — + g 
in Eqs. (O) and ( |4.2| ). A short calculation similar to the 
one presented for $0 = leads to the magnon dispersions 



w±(q) 



q2 + i (72 ^ ^^4 + i6A2(e • q) 



-,1/2 



(4.8) 



where e = sin<I>o ei +cos$o ^2 is the unit vector obtained 
by rotating 62 with an angle — <i>o- The emerging picture 
is yet another manifestation of the peculiar nature of the 
C/(l) symmetry (3.4), in some respects similar to the bi- 
section rule discussed in the concluding paragraph of Sec. 
III. In any case, the main conclusion of the present sec- 
tion persists; namely, the acoustical mode develops max- 
imum instability along the direction e and leads to the 
same critical field given earlier in Eq. 
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The nature of the intermediate phase will be discussed 
in Sec. VI. The present section is concluded with a word 
of caution concerning the validity of the continuum ap- 
proximation at nonzero field, which roughly requires that 
gcUsH <^ J. This strong inequality becomes increas- 
ingly marginal for field strengths in the region H > H2- 



V. SPIRAL PHASE 

The calculation of the low-energy magnon spectrum 
in the spiral phase {h < he) is significantly more com- 
plicated, but the general strategy is identical to that fol- 
lowed in Sec. IV. Hence we introduce new fields according 
to 



0.8 



Q = 9 + f, 



$ = 



smf 



(5.1) 



where 6 = d{x) is the profile of the ground-state spiral 
given by Eq. ( p^ ) while / = f{x, y, r) and g = g{x, y, r) 
account for small fluctuations. The special rescaling cho- 
sen in the second equation is equivalent to working in a 
rotating frametB whose third axis is everywhere parallel 
to the direction of the background staggered magnetiza- 
tion n — (sin 0, 0, cos t 



The new fields (5.1) are int roduce d in the complete 
Lagrangian given by Eqs. (4T) and (4^) which is then 
expanded to second order in / and g. The required al- 
gebra is lengthy but the final result for the linearized 
equations is sufficiently simple: 



/-A/ + [/i/ 
g- Ag + U2 g 

where 



2h cos 6 g + 2X sin 6 d2g, 
-2/icos6l /-2Asin6' 82/, 



Ui = - 72cos(26l), 

U2 = 2AV<52 +72cos20 - 272cos26i - 6^ 



(5.2) 



(5.3) 



are effective potentials that can be calculated for any 
desired set of parameters, as explained in Sec. III. The 
general idea that the calculation of the spectrum in a 
spiral antiferromagnet can be reduced to a Schrodinger- 
like problem in a periodic potential is not newllj, but the 
specific structure of Eqs. ( |5.2| ) requires special attention. 

We found it instructive to consider first the special case 
of spin-wave propagation along the x-axis {d2f = = 
d2g) at zero external field {h = 0). This is actually the 
only case for which the low-energy spectrum was previ 
ously calculated starting from the discrete Hamiltonianl 
If we further perform the temporal Fourier transforma- 
tion with frequency lu, Eqs. (5.2) reduce to 



-g" + U2g = iu^g, (5.4) 



where the prime denotes differentiation with respect to 
X. Therefore, in this special case, the eigenvalue problem 
is reduced to two decoupled ID Schrodinger equations of 
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FIG. 6. Magnon spectrum for spin-wave propagation along 
the X-axis at zero field. Solid and dashed lines distinguish 
between acoustical and optical modes, (a) The spectrum in 
a reduced-zone scheme, and (b) the same spectrum in an ex- 
tended-zone scheme including two replicas of the acoustical 
mode centered at ±(". 



the standard type with potentials U\ and U2 calculated 
at zero field. Also note that both potentials are peri- 
odic functions of 26* and thus their period is actually L/2 
where L is the period of the background spiral. 

The eigenvalue problems (|]^) are solved in Appendix 
A. The numerical procedure yields eigenfrequencies uj = 
tj(gi) as functions of Bloch momentum q\. The latter 
can be restricted to the zone [—i'K/L, 27r/L], because the 
period of the potentials is or to the zone [— C, C] in 
relative units defined as in Sec. IV. Several low-lying 
eigenvalues are illustrated in Fig. 6a using a reduced- 
zone scheme. Solid and dashed lines corres pon d to the 
first and second eigenvalue problem in Eq. ( x4 ) and are 
superimposed in the same graph for convenience. We 
also find it convenient to refer to the two types of modes 
as acoustical and optical. In either case, there is only 
one discernable gap that occurs between the first and 
the second band at the zone boundary. The calculated 
boundary gaps are 0.123 meV and 0.049 meV, respec- 
tively, while the absolute gap of the optical mode at the 
zone center is 0.170 meV. All of the above theoretical pre- 
dictions agree with those obtained in Ref. 5 by a different 
method. They also agree with experiment, except for the 
small (0.049 meV) gap that has not yet been resolved at 
zero field. 

The same results are depicted in Fig. 6b using an 
extended-zone scheme. In fact, this figure displays two 
replicas of the acoustical mode centered at ±C. The need 
for two replicas follows from the structure of dynamic cor- 
relation functions in the laboratory frame, rather than in 
the rotating framo^ctually used in the calculation of the 
magnon spectrumtl. Our results in Fig. 6b are obviously 
consistent with both the experimental and the theoreti- 
cal results obtained in the above reference at zero field. 



We are now in a position to extend the calculation to 
the general case of nonzero field and arbitrary direction 
of spin-wave propagation. The external field enters Eqs. 
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FIG. 7. Magnon spectra for spin-wave propagation along 
the a;-axis, at four values of the applied field H , using a highly 
reduced zone scheme. 



( p. 2]) in two distinct ways. First, it affects the structure of 
the potentials Ui and U2 because the background spiral 
is further distorted. Second, the field induces first-order 
time derivativ es w hich originate in the "nonrelativistic" 
term of Eq. (4.1) and couple the two linear equations 
(5.2). Additional coupling between the two equations 
appears in the case of arbitrary direction of propagation 
because 82/ and no longer vanish. Altogether we are 
faced with a nonstandard eigenvalue problem that is also 
solved in Appendix A. 

Here we present explicit results for four typical values 
of the rationalized field h — 0, 0.3, 0.6 and 0.9 which will 
be quoted from now on by their rounded physical val- 
ues H = 0, 0.5, 1 and 1.5 T. In Fig. 7 we illustrate the 
calculated spectrum for spin-wave propagation along the 
X-axis (^2 = 0) using a highly reduced zone scheme. An 
important check of consistency is provided by the fact 
that the H = results of Fig. 7 agree with those pre- 
sented earlier in Fig. 6a, except that the zone is now 
reduced down to [— C/2, C/2] for reasons explained in Ap- 
pendix A. Furthermore, we no longer employ solid and 
dashed lines to distinguish between acoustical and optical 
modes. Such a distinction is not a priori possible in the 
current algorithm because of the coupling (hybridization) 
of the two types of modes at nonzero field. 

One should keep in mind that the extent of the zone 
[— C/2, C/2] slides with the apphed field, a feature that is 
not apparent in Fig. 7 because the scale of the abscissa 
is adjusted accordingly. The incommensurability param- 
eter C — 0.0273 measured at 77 = is used as input in 
our calculation. The calculated values for H — 0.5, 1 and 
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FIG. 8. The low-energy spectra of Fig. 7 using an ex- 
tended-zone scheme as in Fig. 6b. The energy values at the 
five characteristic spectral points denoted by 1, 2, 3, 4 and 5 
are given in Table I. 



1.5 T are C = 0.0271, 0.0264 and 0.0245. 

At first sight, it would seem difficult to extract useful 
information from the highly convoluted spectra shown 
in Fig. 7. Nevertheless, the most vital information con- 
cerning the low-energy dynamics is easily abstracted from 
Fig. 7 because the low-lying bands are clearly segregated. 
In particular, it is still possible to distinguish between 
the acoustical and the optical mode, at least in an oper- 
ational sense. Thus we unfold the first six branches back 
to the zone [— C, C] and then proceed to the extended-zone 
scheme of Fig. 6b including two replicas of the acoustical 
mode centered at ±C. The resulting low-energy spectra 
are shown in Fig. 8. 

The H — entry of Fig. 8 is but a magnified ver- 
sion of the lower-central portion of Fig. 6b, as expected. 
This version reveals a certain "anomaly" that is not con- 
spicuous in Fig. 6b, namely, a relative crossing between 
the two modes in a narrow region around the zone center. 
The calculated maximum splitting of 0.005 meV is within 
the error margin of the continuum approximation and, in 
any case, beyond experimental detection. But the resolu- 
tion of this theoretical curiosity is interesting: when the 
direction of spin-wave propagation departs slightly from 
the X-axis (52 7^ 0) and/or a finite field is turned on, 
the crossing points become avoided crossings. Therefore, 
strictly speaking, the solid and dashed lines must be in- 
terchanged in the narrow region between the two crossing 
points. This explains the apparent slight inconsistency in 
the labeling of the five characteristic points of the spec- 
trum denoted by 1, 2, 3, 4 and 5 in Fig. 8. The calculated 
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TABLE I. Energy in units of meV at the five characteristic 
points of the spectrum denoted by 1, 2, 3, 4 and 5 in Fig. 8. 



H [T] 


El 


E2 


Es 


E4 


E5 


0.0 


0.170 


0.176 


0.266 


0.298 


0.314 


0.5 


0.122 


0.223 


0.276 


0.299 


0.307 


1.0 


0.073 


0.265 


0.286 


0.300 


0.305 


1.5 


0.025 


0.255 


0.285 


0.298 


0.346 



magnon energies at those points are summarized in Table 
I. 

We now concentrate on the optical mode. The gap 
E2 = 0.176 meV calculated at zero field agrees with 
the measured 0.18(1) meV. Our calculation further shows 
that the above gap evolves quickly with increasing field to 
reach the asymptotic value 0.26 meV around which it os- 
cillates mildly. The complete optical mode evolves into a 
snake-like dispersion with energy values in the range 0.25 
meV < E < 0.29 meV. These predictions are generally 
consistent with experimentij. However, some of the finer 
details deserve closer attention. The calculated energy at 
point 5 in the spectrum remains practically constant at 
i?5 ^ 0.31 meV for ^ 1 T, while a steep crossover takes 
place for higher field values which leads to ^ 0.35 meV 
for H = 1.5 T. These predictions are also in agreement 
with experiments. But the calculated splittings of the op- 
tical dispersion E^ — E^ = 0.02 meV and 0.06 meV, for 
H = 1 and 1.5 T, disagree with the measured 0.05 meV 
and 0.11 meV. It appears that the observed splittings are 
better described by E5 - E2 = 0.04 meV and 0.09 meV. 
In fact, the above identification may not be completely 
arbitrary. For instance, the lowest branch in the optical 
dispersion measured for _ff = 1.5 T shows a clear local 
maximum of 0.28 meV at the zone center, which agrees 
with the calculated maximum E3 = 0.285 meV at the 
zone boundaries ±C rather than the gap E2 = 0.255 meV 
at the zone center. It seems that the lowest branch in the 
observed optical dispersion for = 1.5 T is composed 
of two replicas of the calculated dispersion gentered at 
±C. On the other hand, experimental dataQ at higher 
energies not shown in Fig. 8 indicate the appearance of 
two replicas centered at ±2^. Unfortunately, we cannot 
resolve this issue of proper replication of the basic modes 
because our current formalism does not directly address 
the relevant dynamic correlation functions. 

Next we discuss the acoustical mode. Our calculation 
shows that the energy at point 4 in the spectra of Fig. 8 
remains remarkably stable at E4 w 0.30 meV for all field 
values considered. This feature is also in agreement with 
experiment which indicates only a mild decline from the 
above value with increasing field. Nevertheless, a clear 
disagreement occurs in the lowest branch of the acoustical 
mode. Although explicit cLata points are not given for this 
branch by Zheludev et al.U, the solid lines in their Figs. 
6 and 7, and the corresponding wording in their text, 
suggest that the lowest branch in the measured spectrum 
is also largely insensitive to the applied field. In contrast. 



our calculation predicts a robust reduction of the energy 
gap El with increasing field (see Fig. 8 and Table I). The 
calculated spin- wave velocity is also reduced, albeit at a 
slower rate. 

The preceding apparent disagreement with experiment 
is especially important because it is directly related to 
the issue of local stability of the spiral phase. Indeed, a 
careful numerical investigation reveals that the gap Ei 
vanishes at the critical field hi « 1.01, or Hi w 1.70 T, 
while an unstable mode develops for H > Hi. This mode 
is first detected by the appearance of a real eigenvalue in 
the matrix M of Eq. (A5), when H crosses Hi, which 
corresponds to purely imaginary frequency. As the field 
increases beyond Hi the instability occurs over a non- 
trivial region in q-space. Therefore, the flat spin spiral 
constructed in Sec. Ill is predicted to be locally stable 
only for H < Hi < He- rnp 

It is interesting that the experimental workclu already 
provided evidence for the existence of a critical field 
Hi = 1.7 T that coincides with our theoretical predic- 
tion. However, one should also contemplate the possi- 
bility that such a coincidence may be fortuitous, in view 
of the apparent contradiction between experimental and 
theoretical predictions for the gap Ei. In any case, our 
current result together with the discussion of Sec. IV 
clearly suggest the existence of an intermediate phase in 
the field region 1.7 T < < 2.9 T. The nature of the 
intermediate phase is discussed in Sec. VI. 

In the remainder of this section we take a different 
view of the low-energy magnon spectrum by considering 
spin-wave propagation along the normal to the plane of 
the flat spiral. Our algorithm is adapted to this case 
simply by setting the Bloch wave number qi = and 
calculating frequencies as functions of the wave number 
52 in the j/-direction. It is interesting that no theoretical 
or experimental results exist in this case even at zero 
field. Our results are illustrated in Fig. 9 for the same 
set of field values employed in the preceding discussion. 

The most stable feature of Fig. 9 is its lowest branch 
which exhibits quadratic dependence on q2 near the ori- 
gin. Clearly this branch is the extension of the acousti- 
cal dispersion in the y-direction originating at its points 
where E = 0. Therefore, the complete acoustical mode 
is Goldstone-like in the x-direction but ferromagnetic-like 
in the y-direction. Such a characteristic anisotropy is in 
some respects similar to the situation encountered in the 
spin-flop phase discussed in Sec. IV. 

Higher branches labeled as 1, 2, 3, 4 and 5 in Fig. 9 also 
possess a simple interpretation, for they are the exten- 
sions in the y-dircction of the special spectral points num- 
bered accordingly in our earlier Fig. 8. In contrast to the 
fundamental ferromagnetic-like branch, higher branches 
evolve vigorously with the applied field. In particular, 
branch 1 in Fig. 9 is quickly depressed with increasing 
fleld to become degenerate with the fundamental branch 
at the critical field Hi = 1.70 T not included in the fig- 
ure. For H > Hi this mode becomes unstable over a 
nontrivial region of wave numbers around the origin. Of 
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FIG. 9. Magnon spectra for the same field values as in 
Figs. 7 and 8 but spin-wave propagation along the j/-axis. 



course, this is the instability described earlier in the text 
viewed from a different perspective. 

We have thus provided a fairly complete theoretical 
picture of the low-energy magnon spectrum, including 
predictions for which there exist no experimental data 
at present. It is interesting to see whether or not fu- 
ture experiments could resolve the apparent discrepancy 
in the field dependence of the magnon gap Ei , and thus 
illuminate the important issue of local stability of the spi- 
ral phase, as well as confirm the predicted characteristic 
anisotropy in the low-energy spectrum. 



VI. INTERMEDIATE PHASE 

We now focus on the predicted intermediate phase and 
examine its nature through a direct numerical mini miza - 
tion of the complete energy functional W of Eqs. (3.1) 
and (3.2). The method of calculation is a relaxation al- 
gorithm formulated on the basis of a discretized form of 
the energy functional defined on a square grid. After 
long experimentation with 2D simulations, it progres- 
sively became apparent that the optimal configuration 
for h > hi is actually a ID nonflat spiral characterized 
by a staggered magnetization whose three components 
are all different than zero. 

Therefore, an accurate calculation of the nonfiat spi- 
ral was eventually obtained by a relaxation algorithm 
applied directly to a ID restriction of the energy func- 
tional whose variation leads to the coupled stationary 
equations: 



a> 

CO 

E 

T3 
0) 

O 
O) 
D) 




FIG. 10. Profile of the nonfiat spiral for h = 1.21. The 
three curves correspond to the three components of the stag- 
gered magnetization ni, n2 and 713. The calculated period is 
L = 8.84. 

e" -I- (72 - $'2) cos e sin e -2 A sin^ 6 sin $ 

(sin^e^')' = 2A sin^e sin$ e' (6.1) 

These are ordinary differential equations because both 
angular variables 8 and <i> are assumed to be functions of 
the single coordinate x, while the prime again denotes dif- 
ferentiation with respect to x. Nevertheless, it does not 
seem possible to obtain analytical solutions of Eqs. (|6T|), 
except for the case of the flat spiral ($ = 0) discussed in 
Sec. III. A significant obstacle is the fact that the period 
of the nonflat spiral is not known a priori. Hence our 
numerical solution was carried out on a periodic ID grid 
with specified length L, until a relaxed configuration was 
obtained with energy density w — w{L). We then varied 
L to achieve the least possible energy for each field h, 
and the corresponding optimal period L = L{h). 

An important check of consistency is that the above al- 
gorithm reproduces the results for the flat spiral obtained 
more directly in Sec. Ill, but only when h < hi — 1.01. 
Instead, a nonflat spiral emerges as the optimal solution 
for h > hi. The calculated conflguration is illustrated 
in Fig. 10 for a field value h = 1.21 deliberately chosen 
to be equal to the critical field he of the conventional CI 
transition. The energy of the nonflat spiral depicted by a 
dashed line in Fig. 3a is smaller than the energy of both 
the flat spiral and the uniform spin-flop state throughout 
the intermediate region hi < h < h2- One should also 
stress that the nonflat spiral is here predicted to occur 
for a field applied strictlw along the c-axis, and is not due 
to sample misaljgnmentu or the presence of a transverse 
magnetic fields. 

In a sense, the predicted intermediate phase smooths 
out the original sharp CI transition. This smoothing is 
also apparent in the calculated field dependence of the pe- 
riod L = L{h) which is inserted in Eq. (3.17) to yield the 
results for the incommensurability parameter shown by a 
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dashed line in Fig. 3b. The same figure displays experi- 
mental data taken from Ref. 4 where they were analyzed 
in terms of the conventional CI transition based solely on 
a flat spiral. It should be noted that both the measured 
zero-field incommensurability parameter <^(0) — 0.0273 
and an experimental critical field = 2.15 T were used 
as adjustable parameters in the theoretical analysis of 
Refs. 4 and 5 to obtain a reasonable overall fit. Yet 
the experimental data indicate some smoothing of the CI 
transition near the critical field. This fact is made ap- 
parent in our Fig. 3b where theoretical results for both 
the flat spiral (solid line) and the nonflat spiral (dashed 
line) are calculated using as input only the zero-field pa- 
rameters given earlier in Eq. (3.16). 

Nevertheless, the results of Fig. 3b cannot be inter- 
preted as unambiguous evidence for the existence of an 
intermediate phase, especially because the experimen- 
tal data were taken at the relatively high temperature 
T = 2.4 K. It is feasible that the T = theoretical pre- 
dictions could be further focused by invoking deviation 
from KSEA anisotropy that is allowed by symmetry; i.e., 
by repeating the calculation for nonzero values of the free 
parameter k. One should also keep in mind that a com- 
pletely accurate description of the CI transition may not 
be attainable within the classical approximation. 

The nonfiat spiral exists as a stationary point of the 
energy functional throughout the intermediate phase and 
degenerates into a uniform spin-flop state polarized along 
the y-axis near the upper critical fleld hi = 1.73. Actu- 
ally, our calculation was not pushed all the way to the 
critical field hi because of numerical difficulties that oc- 
cur as the period grows to infinity. The theoretical anal- 
ysis should be completed with a detailed study of the 
stability and dynamics of the nonfiat spiral within the 
full 2D context, in a manner analogous to our treatment 
of the flat spiral in Sec. V. The required computational 
effort is too great to be included in the present paper, 
especially because the proflle of the nonflat spiral is ob- 
tained numerically through the relaxation algorithm. A 
future analysis could, in principle, reveal the existence of 
yet another critical fleld within the intermediate region, 
beyond which the nonflat spiral may cease to be locally 
stable. It is thus important to also examine the nature 
of instability at the upper critical fleld hi, as discussed 
further in Appendix B. 

The conflguration of Fig. 10 may be viewed as a coni- 
cal spiral that nutates around the y-axis. It is interesting 
that a simple conical spiral without nutation had been 
discussed theoretically in connection |jjpt|h the cholcstcric- 
nematic transition in liquid crystalsOtiJ but has not yet 
been observed experimentally because its realization re- 
quires an anomalously small bend moduluslliJ. In con- 
trast, the parameters of Ba2CuGe207 favor the occur- 
rence of the currently predicted intermediate phase. 



VII. CONCLUSION 

We have presented a fleld theoretical description of 
the low-energy dynamics in the spiral antiferromagnet 
Ba2CuGe207. We have thus been able to calculate the 
low-energy magnon spectrum for any strength of the ap- 
plied fleld and any direction of spin-wave propagation. In 
this respect, the present work signiflcantly extends the re- 
sults of Ref. 5 where the spectrum was calculated only at 
zero field and for propagation along the direction of the 
spiral. Therefore, our theoretical results are relevant for 
the analysis of experimental data obtained for nonzero 
field, which were previously analyzed mostly in terms of 
empirical formulas. 

An interesting byproduct of this detailed spin-wave 
analysis is the identification of the two new critical fields 
Hi and Hi, and a corresponding prediction of an inter- 
mediate phase that does not seem to be inconsistent with 
available experimental data. The apparent discrepancy 
in the field dependence of the magnon gap Ei pointed 
out in Sec. V needs to be clarified, but could be due 
to poor experimental resolution at this rather low en- 
ergy scale (0.1 meV or less). The field dependence of 
the incommensurability parameter discussed in Sec. VI 
could be rectified by invoking a slight deviation from the 
KSEIA limit that is allowed by symmetry. Susceptibility 
dataci taken at T = 2 K display a rounded maximum 
which could be explained as a finite-temperature effect 
but does not a priori exclude an intermediate phase. Fur- 
thermore, the set of data for the magnon dispersion dis- 
cussed in connection with Fig. 5 is too limited to provide 
a clear picture. Therefore, a clear identification or dis- 
proof of the intermediate phase may require additional 
experimental work guided by the theoretical predictions 
of the present paper. 

On the other hand, it is desirable to carry out a com- 
plete theoretical analysis of the stability and dynamics 
of the intermediate phase along the lines outlined in Sec. 
VI. A related project is to extend our approach to the 
case Cj£ a field applied in a direction perpendicular to the 
c-axiscl. The field-dependent modifications of the spiral 
can be computed on the basis of Eq. ( t?.19| ) , and a corre- 
sponding calculation of the low-energy magnon spectrum 
can be carried out by a straightforward extension of the 
methods developed in Sec. V. 

Finally, we must comment on the two basic approxi- 
mations made in the present work. The adopted classical 
approach is equivalent to the usual semiclassical approxi- 
mation obtained by the 1/s expansion restricted to lead- 
ing order. The omitted quantum (anharmonic) correc- 
tions are not negligible in this 2D problem but are offset 
in part by the fact that the input parameters are consis.- 
tently estimated within the classical approximationETEl. 
One should also question the validity of the continuum 
approximation whose relative accuracy can be roughly 
estimated from « 0.03 at zero field, but may deterio- 
rate in the presence of a strong external magnetic field. 
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Incidentally, the corresponding parameter e in a typical 
weak ferromagnet such as an orthoferrite (YFeOa) or a 
high-Tc superconductor (La2Cu04) is at least one order 
of magnitude smaller. In any case, the physical picture 
derived is sufficiently complete to provide a basis for a 
meaningful discussion of further refinements. 
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APPENDIX A: EIGENVALUE PROBLEMS 



The eigenvalue problems (5.4) were solved numerically, 
as explained here for the first equation. Taking into ac- 
count that the period of the potential is the Bloch 
representation of the wave function reads 



(Al) 



and the wave equation becomes 

oo 

fm = C^Vn, (A2) 



m— — OO 



where the Fourier coefficients of the potential are given 

by 



C/i. 



L/4 



'-L/4 

4 

— / ' 



exp{—iAn'Kx I L) Ui[6(x)\ dx 
Ui{e) dd 



Ami , , 



(A3) 



Here we use the fact that l]\ is an even function of Q or 
X, and X — x{9) is g iven by the integral ( |3.8| ). Thus the 
last step of Eq. (A3) is in effect a double integral that is 
computed by an ada ptiv e Newton-Cotes algorithm. The 
eigenvalue equation (A2) is then solved by diagonalizing 
the finite matrix that results from a restriction of the 
indices m and n to the interval [—N, N] where N can 
be as low as 20. To be sure, only the first few Fourier 
coefficients of the potential Ui and U2 are important, as 
demonstrated in Table II using as input the zero-field 
parameters quoted in Sec. III. The numerical procedure 
just described yields eigenfrequencies lo = w{qi) as func- 
tions of Bloch momentum qi that can be restricted to the 



TABLE IL Fourier coefficients of the potentials Ui and 
U2 at zero field. The table should be completed with the 
symmetry relations Ui,-n = Ui,n and 172, -n = U2.n- 



n 


Ui,„ 


C/2,„ 





0.13034455 


0.53189772 


1 


-0.49358342 


-0.24049378 


2 


-0.06461030 


-0.04790238 


3 


-0.00637043 


-0.00527215 


4 


-0.00055833 


-0.00048614 


5 


-0.00004588 


-0.00004113 


6 


-0.00000362 


-0.00000331 


7 


-0.00000028 


-0.00000026 


8 


-0.00000002 


-0.00000002 



zone [— 27r/L, 27r/L], or to C] i^i relative units defined 
as in Sec. IV. 

We now return to the general case of nonzero field and 
arbitrary direction of spin-wave propagation. We first 
rewrite Eqs. (5.2) in a form that contains only first- 
order time derivatives. Hence we treat u = f and v = g 
as independent fields and introduce the four-component 
spmor X defined from X'^ = (u, w, /, g). Then Eqs. (O) 
read 



X = MX, 
where M is the differential operator 



M 



/ 

/ 

-Di D3 D4 

-Ds -D2 -Di 



(A4) 



(A5) 



Here Di = -A + C/i, = -A -h U2, D3 = 2Xsm0 62, 
D4 — 2/1 cos 6*, and / is the unit operator. The chief ad- 
vantage of M is that it does not contain time derivatives. 
A superficial disadvant age is that M is not a hermitian 
operator. In fact, Eq. (A4) suggests that the eigenvalues 
of M are purely imaginary and come in pairs ±iLU where 
uj is the desired physical frequency. A real eigenvalue in 
M would correspond to purely imaginary physical fre- 
quency and thus indicate instability of the ground-state 
spiral. All of these features are explicitly realized in the 
following numerical calculation. 

Our task is then to construct a matrix representation of 
the differential operator M. Attention should be paid to 
the fact that the Bloch theorem must now be applied with 
the f ull p eriod L of the spiral because of those terms in 
Eq. (A5) that arc proportional to cos 61 and sin^. Hence 
the operator Di ~ —A + l]\ is replaced by a matrix 
{D\,nm) with elements 



i*!,™™ = (gi + 2n7r/L)^ + 92 



(A6) 



where q\ is now restricted to the zone [— tt/L, tt/L], or 
[— C/2,^/2] in relative units, while 52 is unrestricted be- 
cause the spiral depends only on x. Accordingly, the 
Fourier coefficients of the potential are given by 
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COS 



2mT 



i9) 



Ui{9) d0 



7^008^6 



(A7) 



which differs from Eq. (A3) only in that the fuU period 
L, instead of L/2, is employed. As a result, odd coef- 
ficients in Eq. (A7) vanish, while the collection of even 
coefficients coincides with that obtained from Eq. (AS). 
The operator D2 is treated in exactly the same way re- 
placing Ui with U2- On the other hand, the operator 
= 2Asin6' in Eq. ( [A5| ) is replaced by ^Xq^S where 
S is an antisymmetric matrix whose n-th codiagonal has 
all its elements equal to 



x{0) 



sin 9 dO 



\J 6"^ + 7^cos-^6' 



(A8) 



and Z?4 — 2hcos9 is replaced by 2hC where C is a sym- 
metric matrix whose n-th codiagonal has all its elements 
equal to 

^"■T'' cos^ d9 

—x{9) 



C 



COS 



■\/P + 7^cos^ 



(A9) 



An interesting fact is that both Sn and C„ vanish for even 
n. The most important terms are those with n = ±1, 
whereas higher-order terms account for distortion of the 
spiral from its ideal shape 9 = Ax. Such a distor- 
tion occurs even at zero field in the presence of KSEA 
anisotropy. 

A finite-matrix representation of the differential op- 
erator M is then obtained by restricting the indices m 
and n to the finite interval [—N^N] where iV may again 
be chosen as low as 20. The resulting nonsymmetric 
4(2iV -I- 1) X 4(2A^ -|- 1) matrix is diagonalized numeri- 
cally to yield eigenvalues that are indeed purely imagi- 
nary and come in pairs zbiw where w = w (91,(72) is the 
sought after physical frequency. We have thus obtained 
a number of results using as input the spiral parameters 
A = 1, 7^ = H- /i^, (5 = 5{h) and L = L{h) calculated for 
each field h as explained in Sec. III. The numerical bur- 
den is insignificant and can be carried out interactively. 
Explicit results are discussed in Section V. 



APPENDIX B: VORTEX STATES 

In the original picture of the CI transition! the high- 
field commensurate phase is rendered unstable through 
domain-wall nucleation at the critical field to become 
a spiral phase for h < he- The instability at the higher 
field /i2 > he suggested by the spin-wave analysis of Sec. 
IV is clearly caused by 2D fluctuations. Therefore, it is 
conceivable that the uniform spin-flop phase is actually 
destabilized by nucleation of 2D vortices rather tban ID 
domain walls, as advocated by Bogdanov et a\s3 in a 
number of related models. 

We thus search for genuinely 2D stationary points of 
the static energy that are compatible with U{1) symme- 
try. First, we introduce the usual polar coordinates (r, ij:) 
from 




FIG. 11. The vortex profile 713 — cos6 for three values of 
the parameter 1/ — 0,1 and 2, including the critical value 



X — — cosip, 
7 



— simp, 

7 



(Bl) 



where the overall rescalling by the constant 7 will simplify 
subsequent calculations. A configuration that is strictly 
invariant under the U{1) transformation (3.4) reads 



e^9{r). 



(B2) 



where the minus sign in the second equation is again 
due to the peculiar nature of U{1) symmetry in the 
present problem. Under normal circumstanceSj-e.g., an 
isotropic antiferromagnet in an external fieldllj, both 
choices ^ = tp and ^ = —ip are compatible with axial 
symmetry and are referred to as vortex and antivortex. 
Here only antivortices are possible within the axially- 
symmetric ansatz but wil l be called vortices for brevity. 

When the Ansatz (B2) is introduced in the potential 
V of Eq. (4.2) the corresponding total energy W = 
J Vdxdy reads 



W 



rdr 



sm'^9 



cos 



d9 cos9 sin^A 
dr f ) 



(B3) 



where v = 2 A/7 is the only relevant parameter in this 
static calculation. Also note that we have drqped the 
additive constant term A^/2 from the potential (4.2) and 
thus the energy of the uniform spin-flop state is set equal 
to zero. Variation of the energy functional (B3) with 
respect to the unknown amplitude 9{r) leads to the or- 
dinary differential equation 



dr'^ 



d£ 
dr 



cos9 sin6' — v sin 



(B4) 
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which reduces to the famiHar equation for ordinary spin 
vortices in the extreme limit 1^ = 0. For v ^ 0, solutions 



of Eq. ( B4 ) exhibit slow decay at large distances, namely 



TT 

2 



(B5) 



which turns into exponential decay for = 0. Explicit 
solutions were obtained by a straightforward relaxation 
algorithm and are illustrated in Fig. 11 for three charac- 
teristic values of the parameter 1^= 0, 1 and 2. 



One may restrict the integral in Eq. (B3) to the finite 
range < r < R and examine its behavior for large R. 
A short ca lcul ation taking into account the asymptotic 
expansion (JBq) leads to 



W = n (1 — In R + finite terms, 



(B6) 



and thus the energy exhibits the familiar logarithmic di- 
vergence. This asymptotic result demonstrates the cru- 
cial role played by the parameter u. For u < 1, the energy 
of a single vortex is greater than the energy of the uni- 
form spin-flop state by a logarithmically divergent quan- 
tity. This is the usual situation encountered in the case of 
ordinary vortices {v = 0). The vortex energy is finite for 
1/ = 1 and becomes again logarithmically divergent but 
negative for v > 1. The special point — 2X/j = I leads 
to the same critical field given earlier in Eq. (4.6). 

Therefore, for h < h2, the energy of the uniform spin- 
flop state can be lowered by vortex nucleation. Because 
of the logarithmic dependence of the energy on the size 
of the system, it is clear that a single vortex cannot by 
itself produce a thermodynamically significant effect. In- 
stead, one should expect that a large number of vortices 
is created for h < h2, probably in the form of a vortex 
latticeEil. We have actually performed several numerical 
experiments using the full 2D relaxation algorithm de- 
scribed in the beginning of Sec. VI. Although we have 
already obtained some "spectacular" pictures indicating 
the formation of a vortex lattice, we have not yet been 
able to lower its energy below that of the nonflat spiral. 
It appears that the complete (2D) energy functional dis- 
plays glassy behavior in the intermediate region, which 
may lead to several nearly degenerate local minima. 
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